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Abstract. 

Exact solutions of Schrddinger and Pauli equations for charged particles in an ex- 
ternal stationary electromagnetic field of an arbitrary configuration are constructed. 
Green functions of scalar and spinor particles are calculated in this field. The corre- 
sponding equations for complex energy of particles bounded by short range potential are 
deduced. Boundary condition typical for 5 - potential is not used in the treatment. Ex- 
plicit analytical expressions are given for the shift and width of a quasistationary level 
for different configurations of the external field. The critical value of electric field in 
which the idea of quasistationary level becomes meaningless is calculated. It is shown 
that the common view on the stabilizing role of magnetic field concerns only scalar 
particles. 

I. Introduction. Formulation of the problem. 

The influence of external electromagnetic fields on nonrelativistic reactions with 
charged particles and on the behavior of bound systems (like atoms, ions and atomic 
nucleuses) has being investigated systematically for a long time see, (for example, 
[1-15] ). However, a set of problems, in our opinion, requires additional studying. First, 
it is the appearance of physically meaningless divergences in equations for the complex 
energy of quasistationary bound states. Second, the question about the behavior of 
fermionic systems in super strong electromagnetic fields is not enough elucidated. Note 
that the basic results in this field of physics were received by nonperturbative methods. 
Usually the exact solutions of Schrddinger equations with Hamiltonians taking into 
account the presence of external fields are used [16,17]. It is important, that the 
field can become intensive not only due to increasing of the strength, but also due 
to decreasing of binding energy of the system. In other words, in the vicinity of the 
reaction threshold an arbitrary external field is strong. Third, the common view on the 
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stabilizing role of magnetic field in decays of quasistationary states is inadequate for 
spinor particles [18,19]. That view is caused by the fact that spinor states of electrons 
usually are not taken into accounts in nonrelativistic reactions. In this paper we treat 
a part of these deficiencies. 

Let us consider a charged particle bounded by a short range potential like a 5- 
potential and located in an external stationary electromagnetic field of an arbitrary 
configuration. Note that the potential of zero-radius is a rather widely spread in the 
literature approximation for a multi-electronic atom field and especially for a negative 
ion field, and also for the field of nuclear forces [12,20]. In the general case the 
external field is defined by three independent parameters: strengths of the magnetic H 
and electric E fields and the angle between them. 

The process of particle coming out of a 5-well leads to the transition from an isolated 
energy level to a level in a continuous spectrum, i.e. to a decay of the bound state. 
The external electromagnetic field influences on this transition in two ways. First, it 
forms wave functions of the "free" state and, second, shifts and broads the bound level 
in a 5-well. As the result the initial state of the particle in the external field becomes 
quasistationary. The most adequate instrument for the investigation of such states is 
the complex energy formalism (for example, [4,12,21]) 



where W < is the energy of the nonperturbative bound level, AW is the real shift 
caused by the external field, and T is the width of a level associated with the decay 
probability of the bound state. It is important that a separate consideration of the 
level shift and its width is possible only in a weak external field. In an intensive field 
it is necessary to investigate these quantities in common. 

The main purpose of this article is to derive an equation for the complex energy in a 
stationary external field by the method of an analytical continuation and to solve it for 
some field configurations. This treatment is principally different from the traditional 
derivation using the boundary condition typical for the 5- well [4, 14, 15]. 

The general structure of the paper and its main results are as follows. In the 
second Section explicit solutions for Schrodinger and Pauli equations are constructed 
both for scalar and spinor particles in an external stationary electromagnetic field. 
Then Green functions of the particles are constructed with its help. In spite of rather 
complicated configuration of the external field, these functions can be represented in 
terms of elementary functions. That substantially simplifies their further applications. 

In the third Section the expression for the bound level width in the external field is 
deduced with help of the Green functions. Then the equation for the complex energy 



W = W + AW-iT/2, 
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of the bound particle in the 5-well is built by the method of analytical continuation. 
It has several features different from the ones of the equation traditionally used. First, 
our equation is not transcendental and it is explicitly resolved with respect to the 
complex energy. That is why its right handed side (one fold integral) does not contain 
exponential divergences typical for equations used earlier. Thus it does not require any 
additional regularization. Second, we consistently take into account the spin states, 
which is fundamentally important in the presence of magnetic field. 

Finally in the fourth Section the equation for the complex energy is solved in several 
cases of physical interest. It is well known that the electric field destroys the bound 
state of a charged particle in a potential well. Thus we calculated the critical value of 
the strength of the electric field in which the width of the level becomes comparable 
with its depth in the potential well. In such a field the notion of an isolated level loses 
meaning [21]. When the strengths of the electric and magnetic fields are comparable 
the latter one gives only small corrections to the usual effects: the Stark shift and the 
broadening due to the tunnel effect. But the signs of these corrections are different 
for scalar and spinor particles. In the case when the magnetic field is dominating the 
behaviors of the complex energy of the mentioned particles are principally different. 
For scalar particles the magnetic field exponentially decreases the bound level width 
and therefore stabilizes it. It was assumed that the same effect takes place also for 
electrons (e.g. [5, 10, 14, 15]). But as was shown by the authors of this paper in 
[18, 19], the strong magnetic field linearly increases the level width for particles with 
spin 1/2. Hence in the presence of an external magnetic field it is necessary to take 
into account the spin of the given particle and its contribution to the complex energy 
of the quasistationary level even in non-relativistic problems. 

II. The wave function and the Green function of a charged particle in an 
external stationary electromagnetic field 

In order to receive the wave function of a free charged spinless particle in a stationary 
electromagnetic field of a general configuration, it is necessary to solve the Schrodinger 
equation with the following Hamiltonian: 

h 2 (9 1 

H S h = A _ ihuj H x— — I — moj 2 H x 2 + eEx sin if + eEz cos ip, (2) 

2m oy 2 

where e and m are the charge and the mass of the particle; E and H are the strengths 
of the electric and the magnetic fields, respectively; (p is the angle between the cor- 
responding vectors; ojh = eH/mc is the cyclotron frequency. We suppose that the 
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magnetic field has the orientation along z axis, and the vector of the electric field 
strength lies in the xz plane. 

The Schrodinger equation with the Hamiltonian (2) can be solved by the standard 
method of variables separation. The moving of the particle in a plane normal to 
the magnetic field has the same character as the one in the absence of electric field. 
Corresponding eigenf unctions are Hermit functions, each of them is related to a certain 
Landau level. In this case the coordinate Xq of the orbit centre depends on the conserved 
transverse momentum p y and on the transverse component of the electric field Ecostp. 
The movement of the particle along z axis is similar to the one in the uniform electric 
field, and it can be described by Airy functions [21]. It is not difficult to see that the 
wave function of the charged particle in the field of the considered configuration has 
the form 



h 



(yp y - tw) 



(3) 



ip n (f, t) = Nu n (p) B (f ) exp 
where iV is a normalization coefficient, the Hermit function u n (p) satisfies the equation 

<(p) + (2n + l-pK(p) = 
where n is the number of the Landau level, and the argument of this function is 

(mujH \ 



p(x) 



p v eE sin <p 
x H 7, — 



V h J 

Function B (£) in (3) is proportional to the Airy function Ai regular at infinity 



with the argument 



1/3 



J _. , f'ZmeE cos(p\ , x 

= [ ) \Z-Zo), 

where z is defined by the total energy W of the particle in the external field 

W = hun (n + ^ + eEz cos (p. (4) 

The solution of the Schrodinger equation constructed in that way is valid for a scalar 
particle. In order to take into account spin states of the particle in the concerned non- 
relativistic approximation it is necessary to pass on to the Pauli equation (see, for 
example, [21]) for the spinor wave function ty na (f,t) 
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d - ~ - 

iTl—lpnv (f, t) = H P 1p na (f, t) , 

where the Hamiltonian of the Pauli equation involves the energy of the spin interaction 

— * 

with the external stationary magnetic field H: 

Hp = H Sh + - — aH, 

Zmc 

and a are the Pauli matrixes. Taking into account the concrete orientation of the 
magnetic field this Hamiltonian is reduced to the matrix form 



Hp = 



H Sh + hu H /2 

H Sh - huj H /2 



If we put now into the Pauli equation the spinor wave function of the mixed state 
with an arbitrary time dependence of the coefficients 

where the scalar function ip n (f,t) is defined by the formula (3), we get 



dai/dt = da-ij &t = 0. 

It means, first, that the spin projection on the z axis is conserved. This conservation 
is provided by the homogeneity and the stationarity of the magnetic field. Second, the 
energy of the particle is defined, instead of the expression (4), by the more common 
formula 



W = huuiji + 1/2 + a) + eEz cosy?, (6) 

where a = for the scalar particle and a = ±1/2 for the spinor one. 

Thus four constants (p y , z , n, a) unambiguously parameterize the wave function of 
the particle with spin 1/2 in the external field under consideration. 

Let's construct now a time depending Green function of a charged particle in the 
given field. With the purpose of some simplification of calculations we, without limiting 
the generality, suppose one of arguments of point-to-point Green function to be zero. 
The function will be a diagonal matrix of a kind 
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8 (f, i; 0, 0) = t J / 4*G» (». *. ft) ( r '""exp (i,„ t /2) ) (7) 

n— °-oo -oo ^ ' 

where the scalar Green function Gsh (n,z ,p y ) corresponds to a Schrodinger equation 
G,sh {n,z ,p y ) ~ u n (x)u n (0)Ai(z)Ai(0) exp 

The integral over z included in (7) is easily calculated using the integral represen- 
tation for the Airy function. The integral over p y with Hermit functions is similar to 
the one that was considered in Ref. [22]. The series in Lagerre polynomials formed 
after that is summarized over n with help of the algorithm described in [23] (see also 
[24]). As a result the scalar part of the Green function in a stationary electromagnetic 
field is expressed in terms of elementary functions: 




where 



-foy-tW) . 



mz 



It 



eEzt 



cosy? 



(eEcoscp) t muiH 



1 eE sin tp 
+ -eExt sin ip H 

2 lo h 



24m 

UJ H t 



+ 



cot 



4 

UJ H t 



(x 2 



y ) c °t [~y- 



2xy 



y + 



eEt sin (p 
2mu) 



ii 



(9) 



If we assign in (8), (9) the mass of the particle to unity and perform obvious changes 
in notation we exactly reproduce formulas (A. 7) - (A. 9) from Ref. [14]. At ip = 7r/2 
the relations (8), (9) describe the Green function of a scalar particle in the crossed field 
[9,10,24]. 

3. Derivation of the equation for the complex energy of a bound particle 
in a stationary external field by the method of analytical continuation 

The expression for the probability P of the decay of a particle bound state per unit 
time (or for the width T of the bound level) in the considered field can be received 
with help of the Green function, using the procedure developed in Refs. [11,25]. This 
probability has the form 
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dtG Sh (0, t; 0, 0)/ CT exp (~W t) 

dt f {u H t/2) 



t 1 / 2 sin (u H t/2) 



cxp 



(10) 



where Wo is the energy of the unperturbed bound level (JH), and S at zero coordinates 
is defined by the formula 



(eE cos if) 2 t 3 (eE sirup) 2 1 \uJnt fuJnt\ 



S=~- — + 



24m 2muJ H 



cot 



'Hi 



Poles of the integrand function in points t n = 27m fug are passed from below. 

Function f a (a) is identically equal to unity for scalar particles, but for spinors 
it depends on the polarization of the bound particles in the initial state and of free 
particles in the final one. In the general case 

f a (uiHt/2) = a 1 exp(—iuHt/2) + a\ exp(iuj H t/2), (12) 

where 01(02) are the amplitudes of the probability introduced in (5), that the spin of a 
particle is directed along (or against) the magnetic field. The normalization condition 
demands that a 2 + a 2 = 1. For example, the most realistic situation in the process of 
ionization is the one where bound electrons are not polarized and the detector does 
not distinguish between polarizations of free electrons. In this case for the calculation 
of the total probability of ionization it is necessary to average over polarizations of the 
initial particles and to sum over polarizations of the final ones (see, for example, [21]). 
Then a\ = a 2 . = 1/2 and f a = cos(o>fl-t/2). 

As it is known the width of the bound level is proportional to the probability of 
its decay T = HP. For analytical continuation let us rewrite (10) in such a way that 
the integration in it is carried out on the positive semi-axis. With this purpose it is 
necessary to resolve the exponent in (10) into the real and the imaginary parts and 
in the integral contained cos [(S + W t) /%] to bypass the zero point from below on an 
infinitesimal contour. Such a bypass cancels the divergence proportional to t -1 / 2 of the 
above mentioned integral at zero. As a result the formula for the level width takes the 
form 

1/2 dt 




A\W n \ 1/2 2 \ tt J t 3 / 2 
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2sm(u H t/2) J^. 



a h sin 



S + W t+(-l) k '^t 



+ sin • 



7T 



(13) 



It becomes obvious now that the right hand side of (14) represents the imaginary 
part of the complex expression 



{-W) 



1/2 



i M\ 1/2 7 dt 



2 Vtt 



t 3/2 



2sm(u H t/2) 



— cxp 



Z7T\ 



(14) 



and the left hand side is nothing else than the imaginary part of the expansion (1) 



(-W) 1 ' 2 = (-W - AW + iT/2) 1/2 w | VKq 1 1/2 



AW 



2\W 



1/2 



4|W" 



1/2' 



It is logical to assume that the real part of (14) determines the level shift due to 
the external field. If we use now the identical transformation (recalling that W < 0) 



dt 



exp i 



tWo 

h 



'n\W \ 

h 



1/2 



exp 



/Z7T 

VT 



(15) 



we receive the final equation for the complex energy of a bound level in the stationary 
external field of a general configuration 



-wo 1 / 2 - 



UHt fgjuJHt/2) /iS_ 
2 sm(iu H t/2) 6XP V h 



- 1 



(16) 



where S is defined by formula (11). 

We see that the received equation is resolved explicitly with respect to W. Note, 
that all integrals in intermediate expressions used in the derivation as well as in the 
final formula (16) are finite. Besides the condition of the relative smallness of the 
initial level width and shift has been used. It imposes certain restrictions on the area 
of applicability of (16). However, such restrictions are incorporated in the concept of a 
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quasistationary energy level [21]. More strict mathematical derivation of the equation 
similar to (16) is given in [27]. 

In appendix A of paper [14] a more commonly used variant of derivation of the 
equation for the complex energy in the considered configuration of an external field 
(see also [4,5, 10, 15,20]) is described. Taking into account the spins of particles, we 
can get the equation 



dt ( i TTT 
Wt 



\uut f a (u H t/2) 



2 sm(u H t/2) 



h 

exp(^)-U, (17) 




which differs from the equation (16) only by the substitution Wq — > W in its right hand 
side. Thus the two approaches leading to slightly different equations are formulated 
for the description of quasistationary systems bounded by short-range forces. In this 
connection it is useful to compare the areas of their applicability and to understand 
how much are different their solutions. 

First note that usually an iteration method is used to solve equation (17). Therefore 
the substitution W — > W is accepted on the first step of the iterative procedure, 
reducing thereby (17) to (16). Therefore the received equation (16) is treated usually 
as an approximate one, working only in the case of weak external fields. Thus, the 
choice of one of the alternative approaches suggested above is practically reduced to 
the question: whether iterative corrections of higher orders are necessary or it is needed 
to stop on the "first iteration". 

Obviously, the equation (17) for the complex energy is transcendental and, hence, 
more complicated than (16). Solving a more complicated equation is meaningful only if 
it gives a more exact or an essentially new result. Therefore it is important to discuss 
approximations used by the derivation of the specified relations. At the derivation 
of equation (16) we used conditions T <C |W |, |AW| <C \W \, which impose certain 
restrictions from above on values of external fields. However, the presence of such 
restrictions follows both from the non-relativistic character of the problem, and from 
the initially assumed presence of "the shallow bounded level" in a deep narrow potential 
well [12]. Moreover, the concept of a quasistationary level assumes the smallness of 
its decay probability [21]. We shall discuss the last circumstance in detail below. 

However, it is not possible at all to estimate restrictions on the values of external 
fields, admitted at the derivation of (17). The matter is that all intermediate integrals 
(formulas (A.6), (A.13) - (A.19) in Ref. [14]) as well as the final integral in Eq. (17) 
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are exponentially divergent . At that it is emphasized in works [4,5,20] that the typical 
for a 5-well boundary condition is valid only for a complex value of W 



^w) r% ^-^- ( - w) ^ + (18 > 



Therefore it remains not clear to terms of which order it is necessary to relate the 
last integral at ImW < 0. 

Finally, in our opinion, the key condition 

3/2 

G (f, 0; W) « G (f, 0; W ) « ^ ™ (-W y 1/2 + O (r) . (19) 

on which total procedure of a derivation [14] is based causes some doubt. Namely the 
comparison of asymptotics (18) and (19) gives the traditional equation for a complex 
energy (17). The condition (19) actually demands that not only the singular part but 
also the finite one of the asymptotic expansion of the Green function at zero does not 
depend on external fields. As far as we know, arguments of physical character for the 
benefit of this were not given in the literature. Nevertheless the asymptotic of the 
field Green function at zero (18) contains in the finite part the integral, divergent for 
complex values of W. In these conditions the requirement of the equality of finite terms 
is an obvious excess of the accuracy of the used approximation. This will be shown in 
the following section for the case of a pure electric external field. 

The mentioned above difficulties of the derivation of expression (17) are charac- 
teristic not only for an approach suggested in Ref. [14]. They are typical for any 
derivation of equation (18) using the boundary condition of the (^-potential with a 
complex energy. There are no such problems in Refs. [4,5,10] for the reason that their 
authors determine the complex energy, in contrast to (|T|). with a positive imaginary 
part. It removes the problem of divergences in formulas like (17)-(19), but contradicts 
to the classical definition of the complex energy. In agreement with such a definition 
all probabilities in a quasistationary state should exponentially increase in time [21]. 

Thus, it does not seem to be possible, in our opinion, to derive equation (17) with 
the complex energy W in the right hand side consistently and accurately from the 
mathematical point of view. It is especially problematic to trace all approximations 
used in such a derivation. 
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4. Solving equations for the complex energy of a bound particle for some 

configurations of external fields. 



Let us stop now on some variants of the solution of the received equations for simple 
configurations of an external field. It is stated in Ref. [28] that the elimination of 
divergences in integrals like (18) by the substitution W — > Wo is incorrect. Let's discuss 
further a rather simple configuration of an external electrostatic field from the point 
of view of the validity and the applicability of equations (16) and (17). We use the 
notation of Ref. [28]: the energy e and the strength of the electrostatic field F are 
measured in units of \W \ and in characteristic atom units, respectively: 

W „ eEh 

e = 



i TT »i' ' " H^oi 3 ) 172 



Of course, the particle spin does not show itself in an electrostatic external field. 
At that our equation (16) takes the form [27] 



,1/2 



-e ) 1/2 + 



TT 



At (pi + iAi') + e F- 2/3 Ai (Bi + iAi) 



(20) 



where the common argument of the Airy functions Ai, Bi and of their derivatives 
Ai', Bi' is — SoF~ 2 / 3 , and the zero approximation gives e = — 1. Formula (17), as it is 
affirmed in [15,28], gives in such an external field the transcendental equation 



1 + -F 1 ' 3 \Ai (Bi + iAi') + eF~ 2/3 Ai {Bi + iAi) 

7T L V / 



0. 



(21) 



and the complex energy is contained in arguments of the Airy functions — eF~ 2 / 3 . 
However, note once again that it is impossible to pass from (17) to (21) by means of 
identical transformations. The Airy functions on the left handed side of equation (21) 
are finite while the integral on the right hand side of (17) diverges. It is easy to see, at 
what stage of transformations the infinity "can disappear" . For example, it is possible 
to pass from the general expression (17) to (21) with help of transformation (15) with 
real Wq. However, the integral on the left hand side of (15) diverges explicitly for 
complex W with a negative imaginary part, but the right handed side remains finite. 
The mentioned in the previous section excess of the accuracy in the derivation of 
equation (17) from the comparison of asymptotics (18) and (19) reveals here. 

At the same time, it is possible to receive the transcendental Eq. (21) by analytical 
continuation of the right handed side of (20) over the argument of the Airy functions 
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into the complex plane, replacing e by complex e. As the calculations show, values 
of e determined from (21) are not differ too much from the values that gives the 
explicit formula (20). In a weak field (F< 1) both relations reproduce the well-known 
asymptotic expression [4] 

F 2 i / 4 \ 

if to be limited to first non-vanishing terms of the Airy functions expansions. In the 
strong field F — 1 (i.e. comparable with intra-atomic one) it is easy to get from (20) 

-e « 1,0411 + 0,0451 -i 
and the numerical solution of (21) gives [15,28] 

-e « 1, 0442 + 0, 0388 ■%. 

Besides, there is a critical value of the electric field F cr at which Re(—e) = Im(—e), 
i.e. the energy gap between the "shifted and broadened" level and the continuous 
spectrum disappears. Obviously, that usual concept of a quasistationary level from [21] 
loses the meaning for such intensive electric fields and requires an additional definition. 
In accordance with (20), that critical value is F cr 13, 26. The transcendental equation 
(21) gives F cr ~ 11,38. 

Finelly note that there is a well-known analogy between the complex energy in 
a theory of quasistationary system decays and the complex dielectric permeability in 
physics of semiconductors (for example, [7,8,29,30]). Equations for the complex 
dielectric permeability are derived, as a rule, on the basis of well developed dispersive 
methods. 

One more configuration of an external field which we consider is essentially impor- 
tant to show the role of the magnetic field. For simplicity we put the angle ip = and 
average over polarizations of bound particles. It is convenient to measure the magnetic 
field like the electric one in natural atom units (see, for example, [14]): 

eHh 
n = : -. 

mc | W 1 

If h < F <C 1 it is possible to receive [19] from our equation (16) 

, m/2 F 2 h 2 .F ( h \ ( 46 3 / 2 \ 

(-e) 1/2 -1 = — + + + exp . 

V ; 32 48a 8 V 3aF 2 J ^\ 3F J 
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Here the dependent on spin factors a and b take values a = —2, b = 1 + h/2 for 
a scalar particle and values a = b = 1 for a spinor one. Therefore, the corrections to 
the complex energy due to a magnetic field essentially depend on spins of particles. 
The magnetic field decreases the width of a scalar particle bound level and by that 
stabilizes it. On the contrary, the width of the level of a particle with spin 1/2 increases 
in a magnetic field, i.e. the latter strengthens the destabilizing action of an electric 
field. That difference is still more noticeable in a strong magnetic field [18]. This 
field suppresses exponentially the withdrawal of scalar particles from the 5-well, but 
intensifies it linearly for particles with spin: 

Im(-e)^ ~ fcV* exp (-^^) I ~ ^< 2 (^ 3/2 )- 

If we distinguish particle polarizations after ionization it is not difficult to see with 
help of equation (16) and condition (12) that electrons with spins directed along a 
magnetic field (a\ — 1/2, a 2 , = 0) in the initially non-polarized beam behave like scalar 
particles. Thus, the increasing of the level width in a magnetic field is provided only 
by electrons with spins oriented against the field. 

The reason for the phenomenon described above can be seen from equation (6) for 
the total energy of an electron in an external field. Energies of Landau levels enhance 
with the increase of the magnetic field. That complicates tunneling of a particle from 
the bound state. Only one exception is the ground state of an electron with spin 
directed against the magnetic field. The energy of this state does not depend on H at 
all, but its contribution to the total probability of ionization increases with the increase 
of the magnetic field. 

Therefore for an adequate description of the influence of a magnetic field on the 
behavior of fermionic quasistationary systems it is necessary to take into account spin 
states of fermions even if the considered problem is not relativistic. 
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